Critical exponents in stochastic sandpile models 
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Abstract 

We present large scale simulations of a stochastic sandpile model 
in two dimensions. We use moments analysis to evaluate critical expo- 
nents and finite size scaling method to consistently test the obtained 
results. The general picture resulting from our analysis allows us to 
characterize the large scale behavior of the present model with great 
accuracy. 

Sandpile automata 0] are prototypical models to describe avalanche 
transport processes. All these models show a stationary state that after 
a suitable tuning of the driving fields displays a singular response func- 
tion characterized by power law distributed events. These distributions are 
typically bounded by upper cut-offs related to the system size. In anal- 
ogy with critical phenomena, is possible to define a complete set of scaling 
exponents describing the large scale behavior of these models. 

Despite the large conceptual impact and the huge effort devoted to the 
study of sandpile automata in the last ten years, many basic issues, such as 
the precise values of the critical exponents, the identification of universality 
classes and of the upper critical dimension, still lay unresolved. Theoret- 
ically, many approaches || |5) point out that different sandpile models, 
such as the Bak, Tang and Wiesenfeld (BTW) Q and the Manna || models, 
all belong to the same universality class. Theoretical estimates for critical 
exponents have been provided (especially in Euclidean dimension d = 2) 
by means of different methods ^ [5]], and some exact results || can be 
derived from the Abelian structure of the BTW model. Numerical results 



are difficult to interpret, since different methods of analysis typically yield 
different results ||, |9|, [u], [ll], 12, 13 1. This is probably due to intrinsic scaling 



anomalies and finite size effects present in sandpile models. 

Here, we present very large scale numerical simulations of the Manna 
model||, that is the standard example of a sandpile automaton with stochas- 
tic toppling rule. We show that Manna model can be coherently described 
within a finite size scaling (FSS) framework. Critical exponents are eval- 
uated with great accuracy and the results are confirmed by data collapse 
analysis. 

We consider a two-dimensional square lattice of linear size L and as- 
sociate to each site an integer variable Z{ (energy). At each time step an 
energy grain is added on a randomly chosen site {zi — > Zj + 1). When one of 
the sites reaches or exceeds the local threshold z c = 2 a "toppling" occurs:: 
Zi = Zi — 2 and Zj = Zj + 1, where j represents two randomly chosen nearest 
neighbor sites of site i. A toppling can induce nearest-neighbor sites to top- 
ple on their turn and so on, until all the lattice sites are below the critical 
threshold. This process is called an avalanche. Grains are added only when 
all the sites are below the threshold, which corresponds to a fine tuning of 
the external driving field In addition, the model is conservative and 
energy is dissipated only at boundary sites [l], || . 

Avalanches in sandpile models are usually characterized by three vari- 
ables: the number of topplings s, the area a affected by the avalanche, 
and the avalanche duration t. The probability distribution of each of these 
variables can be described as a power law with a cutoff 

p(x) = x - T *g(x/x c ), (i) 

where x = s,a, t. When the system size L goes to infinity the cutoff x c 
diverges as x c ~ L@ x . Under the finite size scaling (FSS) assumption of 
Eq. (H|), the set of exponents {t x ,{3 x } defines the universality class of the 
model. In order to test the above FSS picture and to find an accurate es- 
timate of the various critical exponents, we perform numerical simulations 
of two-dimensional Manna model with open boundary conditions and con- 
servative dynamics. The lattice size ranges from L = 128 to L = 2048, and 
statistical distributions are obtained averaging over 10 nonzero avalanches. 
The direct numerical determination of the exponents t x from the power law 
behavior of the probability distributions contains intrinsic bias due to the 
lower and upper cut-offs. This makes very difficult to get better than a 



10% accuracy. Extrapolations methods have been devised [13], but the es 



timate of their accuracy is rather difficult. For these reasons, we use the 
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Figure 1: (a) Plot of the moments spectrum for the distribution of toppling 
events s.The linear part has slope 2.73.(b) Data collapse analysis for the 
avalanche size distribution. The values used for the critical exponents are 
t s = 1.27 and S = 2.7. 



more powerful moments analysis on the distribution P(x,L), as also sug- 
gested recently by De Menech et al,[|TTf on the BTW model. We define the 



g-moment of x on a lattice of size L as (x q ) l = J x q P{x)dx. If FSS hy- 



pothesis (Eq. ([[])) is valid, at least in the asymptotic limit (x —■ 
transform z = x/L@ x and obtain 

( X q ) L = L/MS+I-^) f z q+r x g^ dz „ L P x {q+l-r x )^ 



oo), we can 



or in general {x q )i ~ L ax ^ q \ The exponents cr x (q) can be obtained as the 
slope of the log- log plot of < x q >l versus L. Using Eq. we obtain 
(x q+1 ) l I (x q ) l ~ LP" or a x (q + 1) — cr x (q) = X , so that the slope of cr x (q) 
as a function of q is the cutoff exponent j3 x = da x (q)/dq. In Fig.s 1(a), 2(a) 
and 3(a), we show the result obtained from the moments analysis of the 
distribution P(s),P(t) and P(a), respectively. In all cases, we get a clear 
linear behavior starting from q ~ 0.7. For smaller q we observe deviations 
from standard FSS, expected because the integral in Eq. @ is dominated 
by the lower cutoff for small q moments. In particular, corrections to scaling 
of the type {x q )l ~ L ax ^ F(L) are important for q < t x — 1 and when q ~ 
t x — 1, logarithmic corrections give rise to effective exponents up to very large 
lattice sizes. By measuring the slope of the linear part of momentum spectra 
a x{<l)i w e obtain the cut-offs exponents j3 s = 2.73±0.02, (3t = 1.50±0.02 and 
Pa = 2.02 ± 0.02. These exponents[14] are in good agreement with previous 



estimates for the manna model |6[ . If FSS is verified, we can compute the 
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Figure 2: (a) Plot of the moments spectrum for the distribution of avalanche 
time duration t.The linear part has slope 1.50. (b) Data collapse analysis for 
the avalanche time duration distribution. The values used for the critical 
exponents are Tt = 1.5 and j3 a = 1.5. 



exponent t x from the scaling relation (2 — T x )f3 x = cr x (l), that should be 
satisfied for enough large sizes. Using the values of (3 X reported in Table I 
and the values obtained for a x (l) we find the exponents t s = 1.27 ± 0.01, 
Tt = 1.50 ± 0.01 and r a = 1.35 ± 0.01. Also in this case the values are in 
agreement with previous extrapolations or direct measurements^, As 
a final consistency test for the FSS framework we have to verify that we get 
data collapse for the distributions P(x.L) by using the exponents obtained 
from the moments analysis. In fact, the FSS scenario states that rescaling 
q x = x/L@ x and P qx = P(x, L)L^ xTx , the data for different L must collapse 
onto universal curves. In Fig.s 1(b), 2(b) and 3(b), we show that very good 
data collapses are obtained for all distributions. 

In conclusion, we have reported extensive numerical simulations of the 
Manna sandpile model in two dimensions ||. We show that contrary to 
other sandpile models, such as the BTW model JjJ, where it is difficult to 
obtain unambiguous results 11, 13], the FSS assumption is satisfied in the 
Manna model. In this way, we obtain the complete set of avalanche critical 
exponents. Work is in progress to obtain comparable results in d = 3, 4. 
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Figure 3: (a) Plot of the moments spectrum for the distribution of avalanche 
area a. The linear part has slope 2.02. (b) Data collapse analysis for the 
avalanche area distribution. The values used for the critical exponents are 
r a = 1.35 and (3 a = 2.0. 
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